FIOOD  CONTROL  SYSTEM  COMPONENT  OPTIMIZATION-HEC-I  CAPABILITY  SEPTEMBER 


TRAINING  DOCUMENT  NO.  9 


FLOOD  CONTROL  SYSTEM  COMPONENT 
OPTIMIZATION -HEC  -1  CAPABILITY 


mi 


CORPS  OF  ENGINEERS 
?  U.  S.  ARMY 


SEPTEMBER  1977 


jL '  h  '  m*** 

NOV  5  1981 


This  document  has  been  approved 
f.i  p  blic  release  and  sale;  its 
di  t  ibution  is  unlimited. 


THE  HYDROLOGIC 

ENGINEERING  CENTER 

/ 

•research 

•training 

•application 

81  11  03  084 


TRAINING  DOCUMENT 


*  *  -jfe. 

•  it*  „  ^  * 


„,L...  «». afT-i'fr*  Ve*.'*1' 


FLOOD  CONTROL  SYSTEM  COMPONENT  OPTIMIZATION 
HEC-1  CAPABILITY 


October  1974 
(Revised  September  1977) 


The  Hydrologic  Engineering  Center 
609  Second  Street,  Suite  I 
Davis,  California  95616 


Accession  Tor 
NTJS  GBIAI 
DTIC  TAS 

Onauvcuaced 

J’ji'.l  if  leal  Ion — 


- - 

pc trlb’-tlon/ 


b:  ilty  0o'>s 
•Av-:  V  -md.'-.r 

...  •  Lr  i-  -  Til 


□  □ 


TABLE  OF  CONTENTS 


Page 


INTRODUCTION  .  1 

BASIC  EXAMPLE  DESCRIPTION  .  2 

HYDROLOGIC  MODEL  .  2 

ECONOMIC  EVALUATION— EXISTING  CONDITIONS  .  4 

FLOOD  CONTROL  MEASURE  OPTIMIZATION  .  5 

SIZING  RESERVOIR  AND  PUMPING  PLANT— UNCONSTRAINED  .  6 

a.  Detention  Storage  .  6 

b.  Pumping  Plant  .  7 

SIZING  RESERVOIR  AND  PUMPING  PLANT- 

HYDROLOGIC  PERFORMANCE  CONSTRAINED  .  12 

SIZING  RESERVOIR,  PUMPING  PLANT  AND  DIVERSION  .  14 

SIZING  LOCAL  PROTECTION  PROJECTS  .  15 

SIZING  RESERVOIR,  PUMPING  PLANT,  DIVERSION, 

AND  UNIFORM  PROTECTION  LOCAL  PROJECTS  .  17 

OBJECTIVE  OF  THE  FLOOD  CONTROL  SYSTEM  COMPONENT 

OPTIMIZATION  UTILIZING  HEC-1  .  18 

REFERENCES  . 19 

APPENDIX  A  -  INPUT  DATA 

EXHIBITS 

1.  Hydrologic  Model  (Existing  Conditions) 


2.  Multiflood,  Multiplan  Model  (Economic  Evaluation  of 
Existing  Conditions) 

3.  Sizing  Reservoir  and  Pumping  Plant  (Unconstrained) 

4.  Sizing  Reservoir  and  Pumping  Plant 
(Hydrologic  Performance  Constrained) 

5.  Sizing  Reservoir,  Pumping  Plant  and  Diversion  (Unconstrained) 

6.  Sizing  Levee  and  Channel  Modification  (Unconstrained) 

7.  Sizing  Reservoir,  Pumping  Plant,  Diversion  and  Uniform 
Protection  Local  Projects  (Unconstrained) 


t 

i 

* 

i 


LIST  OF  FIGURES 

Figure  Wo.  Title  Page 

1  Study  Area  and  Schematic  Representation  ....  3 

2  Adjustment  of  Component  Size  by 

fJewton-Raphson  Convergence  Procedure . 10 

3  Effect  of  Diversion  on  Flood  Hydrograph  .  15 


111 


INTRODUCTION 


HEC-1  has  been  augmented  to  provide  the  capability  of  automatically 
determining  the  sizes  of  flood  control  system  measures  that  result  In 
maximizing  total  system  net  economic  benefits  subject  to  possible  hydro- 
logic  performance  targets.  The  system  flood  control  measures  that  can 
be  automatically  sized  are: 


Detention  storage  reservolr(s) 

.  Pumping  plant(s) 

.  Diversion(s) 

.  Local  protection (s) ,  l.e.,  channel  modification,  levee, 
floodwall 


This  document  presents  detailed  Illustrated  examples  of  facility 
optimization  using  HEC-1.  The  examples  are  designed  to  assist  in  data 
assembly  and  coding,  output  interpretation,  and  study  management. 

Examples  included  are  constructed  in  building  block  sequence  to 
illustrate  the  relationships  between  the  hydrologic,  economic  and  cost 
data  and  demonstrate  selected  capability.  Examples  Illustrated  include: 
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■-€  Hydrologic  Model  for  existing  conditions. 

Economic  evaluation  of  existing  conditions. 

<*-) Optimization  of  Reservoir  and  Pumping  Plant  with  no 
hydrologic  constraints; 

Optimization  of  Reservoir  and  Pumping  Plant  with 
hydrologic  performance  constraints’ 

I*  ''Optimization  of  Reservoir,  Pumping  Plant  and  Diversion 
(unconstrained)' 

^Optimization  of  local  protection  projects*  levee  and 
channel  modification  (unconstrained): 

-^Optimization  of  Reservoir,  Pumping  Plant  and  local 

protection  projects  with  uniform  local  protection  level. 


The  basic  reference  for  HEC-1  is  the  Users  Manual  listed  as  reference m . 
The  Input  data  supplement,  reference  2,  updates  Addendum  6  of  reference  1 
to  include  the  facility  optimization  capability.  Technical  Paper  No.  42, 
reference  3,  describes  the  conceptual  basis  for  the  optimization  problem  and 
explains  the  characteristics  of  the  flood  control  measures  (except  for  the 
local  protection  capability  that  has  recently  been  added)  and  a  field  app¬ 
lication.  Reference  4  summarizes  various  optimization  algorithms  and  also 
Includes  a  list  of  references  pertinent  to  the  subject  matter  presented 
herein.  Reference  5  describes  In  detail  the  methodology  Involved  In  the 
calculation  of  expected  annual  damages. 


BASIC  EXAMPLE  DESCRIPTION 


The  study  area  lies  In  the  flood  plain  of  a  large  river  and  Is 
presently  protected  (to  a  degree)  by  a  major  levee.  The  levee  greatly 
restricts  outflow  from  the  study  area.  Most  of  the  storm  runoff  (within 
the  study  area)  originates  from  the  higher  elevations  (bluff  areas),  and 
most  flooding  occurs  in  the  lower  reaches  of  the  study  streams.  Develop¬ 
ment  in  the  flood  hazard  areas  consists  of  agricultural  crops,  industrial- 
commercial  areas  and  residential  development.  Figure  1  is  a  general  map 
and  schematization  of  the  example  area. 

Proposals  for  protecting  vulnerable  areas  from  potential  flooding 
include  a  detention  storage  reservoir  at  station  10,  channel  modification 
from  station  10  to  30,  levee  from  station  20  to  30,  flow  diversion  (by¬ 
pass)  from  station  20  to  30,  and  a  pumping  facility  with  forebay  ponding 
at  the  basin  outlet,  station  30  (see  Figure  1-a). 


HYDROLOGIC  MODEL 


The  hydrologic  model  for  existing  conditions  is  needed  to  define 
the  base  hydrology  and  provide  a  mechanism  for  evaluating  the  performance 
of  proposed  alternatives.  Care  must  be  taken  in  developing  the  base  model 
to  assure  that  all  feasible  alternatives  can  be  easily  evaluated  and  that 
the  pattern  hydrologic  event  is  reasonably  representative  for  the  area, 
i.e.,  will  not  bias  evaluation  of  alternatives.  Data  required  for  coding 
the  basic  hydrologic  model  is  given  in  reference  1. 

Since  the  primary  objective  of  this  supplement  is  illustration  of 
flood  control  system  component  optimization,  the  hydrologic  model  has  been 
kept  simple  in  that  discharge  hydrographs  of  a  specific  event  are  read  In 
rather  than  computed  from  rainfall -runoff  relations  during  the  optimization. 
(The  hydrographs  were  essentially  computed  in  a  previous  run).  A  hypo¬ 
thetical  event  was  synthesized  that  ranged  in  frequency  from  the  annual 
event  (1.0  exceedence  frequency)  In  the  upper  basin  reaches  to  about  the 
5-year  event  (.2  exceedence  frequency)  in  the  lower  basin.  Channel  routing 
criteria  has  been  developed  for  the  streams  from  multiple  water  surface 
profile  calculations  and  for  the  restricted  outlet  at  station  30  from  the 
geometry  of  the  outflow  culvert  and  local  topography.  Table  1  (Appendix  A) 
contains  a  tabulation  of  the  hydrologic  data  for  existing  conditions. 

Exhibit  1,  page  1  of  2,  Is  a  listing  of  the  HEC-1  input  data  for  the 
hydrologic  model.  The  hydrologic  simulation  of  existing  conditions  Indi¬ 
cates  that  for  the  selected  event,  the  peak  flows  at  stations  10,  20  and 
30  are  5,370  cfs,  5,370  cfs  and  10,154  cfs,  respectively.  The  maximum 
storage  level  achieved  at  station  30  Is  9,557  ac.ft.  (maximum  storaoe  at 
station  30  not  shown  In  computer  orlntout  Included)  and  the  peak  outflow 
is  1,200  cfs. 
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ECONOMIC  EVALUATION-EXISTING  CONDITIONS 


The  economic  evaluation  for  existing  conditions  provides  the  base 
from  which  economic  benefits  of  alternatives  may  be  evaluated.  The  eco¬ 
nomic  evaluation  of  flood  damages  requires  that  flow-damage-frequency 
analysis  be  performed  to  develop  "expected"  (or  average)  annual  damages. 
Reference  3  and  Addendum  3  of  reference  1  discuss  the  general  application 
of  flood  damage  frequency  analysis  to  flood  alternative  evaluations  and 
describe  the  concepts  embodied  in  HEC-1 . 

The  information  required  (in  addition  to  the  hydrologic  model)  is 
flow  (or  storage)  -  damage  relationships  and  exceedence  frequency  relations 
at  the  damage  centers.  Additional  coding  is  required  to  set  up  the  multi¬ 
plan  feature  of  HEC-1  and  establish  the  range  of  floods  needed  to  evaluate 
the  hydrologic  and  economic  effects  of  alternatives. 

The  damages  in  reaches  1030  and  2030  are  mostly  rural  and  result  from 
overflow  from  the  respective  stream  channels.  Damage  surveys  have  devel¬ 
oped  relationships  between  stage  and  damages  for  these  reaches  for  a  number 
of  categories  of  damages.  Water  surface  profile  studies  developed  rating 
curves  for  the  index  stations  as  shown  on  Figure  1  so  that  flow-damage 
functions,  as  required  by  HEC-1,  could  be  developed.  The  damages  at  loca¬ 
tion  305  are  mostly  urban,  conmercial  and  industrial  (and  are  thus  large) 
and  occur  because  of  ponding  behind  the  levee.  In  HEC-1  storage  is  used 
instead  of  stage  to  represent  level  and  thus  a  storage-damage  function  has 
been  developed  at  this  site.  Storage  is  analogous  to  stage  and  the  func¬ 
tion  is  developed  from  the  usual  stage-damage  relationship  and  a  site  stage- 
storage  relationship. 

The  required  exceedence  frequency  relationships  for  stations  10  and 
20  were  based  on  a  partial  duration  series  analysis  because  significant 
damages  occur  from  events  that  occur  more  frequently  than  the  annual  event. 
These  curves  were  developed  from  regional  relationships  developed  In  other 
studies.  The  required  frequency  relationship  for  station  30  is  storage- 
exceedence  frequency.  This  function  was  derived  by  developing  synthetic 
events  that  would  reproduce  the  regional  curves  at  station  10  and  20,  sim¬ 
ulating  the  hydrologic  operation  of  the  system  for  these  events,  and  plot¬ 
ting  the  resulting  peak  storage  levels  for  these  events  versus  their  ex¬ 
ceedence  frequencies.  Table  2  (Appendix  A)  contains  the  economic  and 
frequency  data  for  the  damage  centers. 

The  determination  of  the  range  of  floods  needed  requires  evaluation 
of  the  exceedence  frequency  relations,  base  hydrology  and  damage  relations. 
The  objective  in  developing  the  range  of  floods  (multi-plan  flood  ratios) 
is  to  provide  for  automatic  revision  of  the  exceedence  frequency  relation¬ 
ship  so  that  expected  annual  damages  can  be  computed  for  alternative  pro¬ 
posals.  The  procedure  used  for  automatically  revising  the  frequency  curve 
is  explained  in  Addendum  3  of  reference  1.  To  accomplish  this,  the 
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ratios  should  develop  floods  that  cover  the  range  of  damaging  floods 
at  all  damage  centers;  in  our  example,  the  range  extends  from  the  six  times 
per  year  event  at  damage  centers  1030  and  2030  to  above  the  .005  event  at 
305.  The  ratios  contained  in  Table  2  (Appendix  A)  when  applied  to  the  syn¬ 
thetic  event  of  the  hydrologic  model  adequately  cover  the  range. 

The  multi-plan  coding  has  been  prepared  for  two  plans,  which  is 
necessary  for  the  optimization  examples  following.  The  two  plans  are  both 
for  existing  conditions  which  is  of  course  redundant.  If  the  multi-plan 
capability  were  being  applied  by  itself,  coding  should  be  for  as  many  alter¬ 
natives  as  is  desired  for  study.  Exhibit  2,  pages  1  and  2,  are  a  complete 
listing  of  data  input  with  notations  as  to  revisions  required  from  the  basic 
hydrologic  model  and  additions  for  the  multi-plan  evaluation. 

The  output  for  a  multi-plan  run  includes  complete  hydrologic  simula¬ 
tion  for  existing  conditions  and  the  proposed  plan  of  improvement  (none  for 
example)  for  each  of  the  range  of  runoff  events  (nine  for  the  example)  and 
integration  of  the  damage  relationships.  The  results  indicate  expected 
annual  damages  under  existing  conditions  are  $33,580,  $33,580  and  $1,110,210 
for  damage  reaches  1030,  2030  and  305,  respectively. 

The  economic  output  (printout  for  station  1030  is  page  3  of  Exhibit  2) 
begins  with  a  printout  of  control  codes  and  includes  (1)  a  listing  of  data 
input  (ECONOMIC  DATA  FOR  STATION  1030  PLAN  1)  which  includes  exceedance  fre¬ 
quency  in  events  per  year,  peak  flow  and  damages,  (2)  computation  of  expected 
annual  damages  (FLOOD  DAMAGES  FOR  STATION  1030  PLAN  1)  which  includes  allo¬ 
cation  of  probability  intervals  (PR0B  INT)  to  the  range  of  flood  events  (FLOW) 
and  incremental  computed  damage  contribution  to  expected  annual  damages 
(SUM,  TYPE  1,  etc.)  that  are  based  on  the  product  of  PR0B  INT  and  damage 
associated  with  FLOW,  and  (3)  the  same  information  for  the  alternative  plan. 

If  the  alternative  plan  had  reduced  annual  damages,  then  the  benefits  (AVG 
ANN  BFT)  would  be  positive  and  equal  to  the  difference  between  PLAN  1  and 
PLAN  2. 


FLOOD  CONTROL  MEASURE  OPTIMIZATION 


The  information  required  in  addition  to  the  hydrologic  model  and 
multi-plan  economic  data  for  flood  control  measure  optimization  are  the  per¬ 
formance  parameters  and  cost  relationships  for  the  flood  control  features 
being  considered.  The  mathematical  structure  for  the  optimization  and  the 
search  strategy  are  discussed  in  detail  In  reference  3.  It  should  be  remem¬ 
bered  (or  understood)  that  economic  optimum  is  achieved  when  the  facilities 
are  sized  such  that  the  computed  difference  between  expected  annual  benefits 
and  expected  annual  costs  is  maximized.  The  solution  may  proceed  uncon¬ 
strained  or  it  can  be  constrained  such  that  a  minimum  hydrologic  performance  at 
specified  control  points  must  be  accomplished  simultaneously  with  the  net 
benefit  maximization. 
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The  general  technique  used  is  to  successively  operate  the  multi -pi an 
simulation  in  a  controlled  fashion  while  automatically  adjusting  component 
sizes  toward  optimum. 


SIZING  RESERVOIR  AND  PUMPING  PLANT  -  UNCONSTRAINED 


The  first  optimization  example  will  be  the  determination  of  the  op¬ 
timum  (economic)  sizes  for  a  reservoir  located  at  station  10  and  a  pumping 
plant  to  be  located  at  station  30  that  discharges  through  (or  over)  the 
levee.  There  is  no  minimum  constraining  hydrologic  performance  required. 
Information  must  therefore  be  assembled  and  coded  that  will  describe,  in  a 
general  way,  the  cost  and  performance  of  the  storage  reservoir  and  a  pump¬ 
ing  facility. 

a.  Detention  Storage. —  The  detention  storage  reservoirs  that  may  be 
considered  with  HEC-1  are  those  for  which  it  is  possible  to  define  the 
operating  characteristics  as  unique  functions  of  the  storage  contents  with¬ 
in  the  reservoirs.  A  reservoir  with  an  uncontrolled  outlet  works  exactly 
meets  this  requirement.  To  provide  capability  for  automatic  adjustment  of 
operating  characteristics  (as  is  required  for  automatic  optimization),  a 
reservoir  is  characterized  by  (1)  the  outflow  characteristics  of  a  low  level 
outlet,  which  is  defined  by  the  centerline  elevation  of  the  outlet  and  an 
orifice  equation  of  the  form: 

Q  =  CA  \f2g  (H)EXP . (1) 

where 

C  =  orifice  discharge  coefficient 
A  =  outlet  area 
H  =  head  on  low  level  outlet 
g  =  acceleration  of  gravity 

EXP  =  exponent  dependent  on  tailwater  conditions,  0.5  if  no 
tail water 

and  (2)  the  overflow  characteristics  of  a  spillway  which  is  defined  by  a 
weir  equation  of  the  firm: 


Q  =  C*LH*  3/2 


(2) 


where 

C*  =  weir  discharge  coefficient 
L  =  length  of  spillway 
H*  =  head  on  spillway 
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and  (3)  the  site  storage  characteristics  which  are  defined  by  an  elevation- 
storage  capacity  relationship.  For  an  index  storage  to  be  optimized,  which 
is  the  storage  at  the  elevation  of  the  spillway  crest,  the  above  relation¬ 
ships  are  merged  to  define  the  reservoir's  outflow  as  a  function  of  the 
storage  level  in  the  reservoir  (Modified  Puls  method  of  routing). 

Two  modes  are  possible  for  a  reservoir  optimization.  In  the  usual 
mode  (for  our  example)  a  reservoir  that  can  be  characterized  by  a  low  level 
outlet  and  an  overflow  weir  as  described  above  will  be  automatically  adjus¬ 
ted  in  its  index  storage  capacity,  along  with  all  other  system  components, 
to  achieve  the  minimum  value  of  the  objective  function  (defined  in  reference 
3).  The  alternative  mode,  not  illustrated,  permits  optimization  of  the 
size  of  the  low  level  outlet  assuming  the  reservoir  does  not  spill,  which 
is  appropriate  for  pondage  in  low  lying  areas. 

The  cost  relationships  for  the  reservoir  in  the  usual  mode  consists 
of  a  capital  cost  function  and  an  associated  capital  recovery  factor  for 
converting  the  capital  cost  to  annual  cost,  and  the  annual  cost  of  operation, 
maintenance  and  replacement  expressed  as  a  proportion  of  capital  cost.  The 
capital  cost  function  includes  land  acquisition  and  construction  costs,  in¬ 
terest  during  construction,  etc.,  expressed  as  a  function  of  the  index  stor¬ 
age  size  of  the  reservoir.  The  capital  cost  for  a  specific  reservoir  size 
being  evaluated  during  optimization  is  interpolated  from  this  function  and 
the  equivalent  annual  cost  is  computed  as  the  product  of  the  capital  cost 
and  the  capital  recovery  factor  for  the  appropriate  discount  rate.  The 
annual  cost  of  operation,  maintenance  and  replacement  is  the  product  of  the 
annual  cost  proportion  and  the  interpolated  capital  cost.  The  total  annual 
cost  of  the  reservoir  is  the  sum  of  these  two  costs.  Table  3  (Appendix  A) 
contains  the  data  describing  the  performance  and  cost  of  the  proposed 
reservoir. 


b.  Pumping  Plant  —  A  pumping  facility  removes  volume  from  the  system 
at  a  rate  equal  to  the  pumping  capacity.  The  performance  characteristics 
of  a  pumping  plant  are  defined  by  an  initial  threshold  water  level  at  which 
the  pump  is  activated  and  the  discharge  capacity  of  the  pumping  facility. 

In  this  analysis,  it  is  assumed  that  water  pumped  from  the  system  does  not 
later  appear  at  other  locations  in  the  system.  The  cost  of  a  pumping 
facility  is  computed  from  a  capital  cost  function  and  an  associated  capital 
recovery  factor  for  converting  to  equivalent  annual  cost,  the  annual  oper¬ 
ation,  maintenance  and  replacement  cost  that  is  a  proportion  of  the  capital 
cost,  and  the  annual  power  cost.  The  power  cost  is  adjusted  if  the  volume 
to  be  pumped  changes  as  the  system  components  sizes  are  being  optimized. 

It  can  be  demonstrated  that  no  matter  the  pumping  capacity,  the  power  costs 
would  not  materially  change  if  the  volume  to  be  pumped  does  not  change. 

The  annual  power  costs  are  therefore  adjusted  only  for  water  that  is  removed 
from  the  system  by  diversions  or  other  pumping  facilities.  The  annual  cost 
is  the  sum  of  the  equivalent  annual  cost,  annual  operation  and  maintenance 
cost,  and  annual  power  cost.  Table  4  (Appendix  A)  contains  the  data  descri¬ 
bing  the  performance  and  cost  of  the  proposed  pumping  plant. 
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The  coding  requires  initial  estimates  for  the  facility  sizes 
(starting  values)  and  a  number  of  control  codes  to  indicate  location 
and  type  of  facility  to  be  sized.  The  starting  values  selected  were 
10,000  ac.ft.  and  4,000  cfs  for  the  reservoir  and  pumping  plant,  res¬ 
pectively.  Exhibit  3,  pages  1  and  2,  are  a  listing  of  the  input  data 
for  this  example  including  notations  of  revisions  and  additions  to  the 
data  required  for  the  multi -pi an  evaluation  example. 

Exhibit  3,  pages  3-43,  are  reproductions  of  the  complete  output 
from  the  optimization  run.  The  output  of  an  optimization  run  includes: 

1.  The  derived  optimum  size  for  each  facility  in  the  system 
included  in  the  optimization  (page  43). 

2.  Complete  hydrologic  simulation  of  the  system  with  and 
without  the  optimally  sized  facilities  for  the  range  of 
floods  processed  (nine  for  this  example)  (pages  6  -  42). 

3.  Economic  expected  annual  damage  analysis  with  and  without 
the  optimally  sized  facilities  for  each  damage  center  in 
the  system  (pages  17,  24  and  41). 

4.  Costs  for  the  derived  system  facilities  (pages  11  and  40). 

5.  A  summary  of  system  cost,  performance  and  net  benefits 
(page  42). 

The  derived  optimum  sizes  are  9,119  ac.ft.  for  the  reservoir  and 
2,885  cfs  for  the  pumping  plant  (summary  page  43).  The  total  capital 
cost  is  $7,497,000  and  system  annual  net  benefits  are  $173,000  (benefit 
cost  ratio  of  1.26).  The  derived  values  were  adjusted  from  the  starting 
values  of  10,000  ac.ft.  and  4,000  cfs  which  corresponded  to  a  capital 
cost  of  $8,740,000  and  system  net  benefits  of  158,000  (page  43).  It  is 
necessary,  in  each  case,  to  test  for  possible  local  optima  in  the  search 
procedure.  This  was  accomplished  by  making  a  separate  run  with  starting 
values  of  3,000  ac.ft.  and  500  cfs  respectively.  The  derived  sizes  were 
6,584  ac.ft.  and  2,835  cfs  costing  $6,591,000  and  resulting  in  annual 
system  net  benefits  of  $199,000.  The  results  indicated  that  a  local  op¬ 
timum  did  exist  such  that  additional  runs  were  made  with  different  initial 
values  until  it  could  be  reasonably  concluded  that  the  prop?**  sizes  were 
6,584  ac.ft.  for  the  reservoir  and  2,835  cfs  for  the  pumping  plant. 

The  hydrologic  performance  can  be  characterized  by  the  "degree  of 
protection"  provided,  i.e.,  the  exceedence  frequency  of  the  threshold  of 
damaging  flow.  At  damage  center  1030,  the  zero  damage  exceedence  frequency 
was  reduced  from  about  the  5  times  per  year  event  to  about  the  annual  event 
(deduced  from  page  17  and  the  additional  runs  made).  Note  that  damages  at 
station  1030  are  quite  small  in  relation  to  those  at  305  and  therefore 
probably  had  very  little  Influence  on  the  determination  of  the  optimum  sizes. 
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At  damage  center  305,  the  frequency  of  significant  damages  was  reduced 
from  about  the  3-year  exceedence  Interval  event  to  about  the  10-year 
event,  which  Incidentally  reduced  expected  annual  damages  by  more  than 
half. 

Detailed  study  of  the  output  can  provide  insight  Into  the  optimi¬ 
zation  methodology  as  well  as  the  sensitivity  of  the  system  performance 
to  a  range  of  facility  sizes.  Pages  3  through  6  of  Exhibit  3  contain 
detailed  output  on  the  progress  of  the  optlmazation.  The  variables  for 
optimization  printed  on  page  3  are  defined  below  and  a  review  of  the 
search  procedure  (reference  3)  and  the  corresponding  results  from  the 
output  are  described. 

Variable  Definition 

NC  =  Counter  denoting  stage  in  search  cycle  (1-3) 

M  =  Variable  that  is  being  adjusted  for  this  cycle 
(corresponds  to  fields  on  J2  card  listed  above 
as  SYSTEM  OPTIMIZATION) 

Ml  =  Next  variable  to  be  adjusted  (optimized) 

VAR(M)  =  Current  value  of  variable  M 

VAR(Ml)  =  Current  value  of  variable  Ml 

OBJ  DEV  =  Used  in  connection  with  hydrologic  performance  con¬ 
straint;  described  in  example  in  next  section 

TANCST  =  Total  annual  cost  of  facilities  at  current  values 

ANDMG  =  Total  annual  damage  for  all  damage  centers  for  facilities 
at  current  values 

0  FTN(NC)  =  Objective  function  that  is  being  minimized;  in  this 
example  it  is  the  sum  of  TANCST  and  ANDMG 

Search  Procedure  (see  reference  3) 

(1)  First,  trial  sizes  of  all  system  components  are  nominated  and 
the  entire  system  is  simulated  in  all  of  its  hydrologic,  costs, 
and  economic  detail  to  calculate  the  value  of  the  objective 
function,  which  for  unconstrained  optimization  is  the  sum  of  the 
equivalent  annual  cost  (TANCST)  and  annual  damage  (ANDMG). 
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The  first  value  (NC=1)  of  the  objective  function  is  1018.883 

(2)  Then  the  size  of  one  component  is  decreased  by  a  small 

selected  amount  (1  percent)  and  the  simulation  is  repeated 
for  the  entire  system  to  compute  a  new  value  of  the  objec¬ 
tive  function.  This  is  repeated  again  resulting  in  three 
unique  values  of  the  objective  function  for  small  changes 
in  the  size  of  one  component. 

The  values  of  the  variable  and  objective  function  are 


NC 

VAR(M) 

0  FTN(NC) 

1  f(X0) 

10000 

1018.883 

2  f(X0  -  AX) 

9900 

1018.205 

3  f(X0  -  2aX) 

9800 

1017.645 

(3)  From  these 

three  values,  an  estimate 

is  made  of  the  compone 

size  that  would  result  in  the  minimum  value  of  the  objective 
function.  The  computation  of  the  adjustment  is  shown  in  Figure 
2  and  proceeds  as  follows: 


Figure  2.—  Adjustment  of  Component  Size  by  Newton- Raphson  Convergence  Procedure 
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inwhich  S'  (xo  "  y)  =  [/(X0)  -  f(Xa  -  AX)](AX)-' . (5) 

r  -  y)  -  £/<xo  - 2AX>  ■  2/<x«  -  AX)  +  /<X.)J<AX)'- . (6) 

and  AX  3  incremental  change  In  X;  X  x  size  of  variable  being  optimized; 
XQ  =  present  size  of  component  X;  and  X*  *  projected  "new"  size  for  X. 

The  calculation  for  adjustment  of  VAR(M)  Is  as  follows: 


■  (*--“)][<■  (*--“)]' 


f'  (X0  -  )  8  (1018.883  -  1018.205)/ A  X  -  0.678/aX . (7) 

f"  (XQ  -  -jr  )  -  [  1017.645  -  2(1018.205)  +  101 8. 883]/ AX2  .118/aX2.  .  (8) 
XQ  =  10000;  AX  =  (.01)  (10000)  *  100 

X*  =  10000  -  Q-678/100,  .  =  9380.  (to  closest  10)  .  .  .(9) 

.118/(100)*  * 

(4)  After  adjustment  of  the  size  of  the  system  component,  the 
entire  system  is  simulated  again  In  detail  to  compute  the 
new  value  of  the  objective  function  and,  provided  the  ob¬ 
jective  function  has  decreased,  the  procedure  then  moves 

to  the  second  system  component  whose  scale  is  to  be  optimized. 

The  output  at  this  stage  reads: 

VAR  1  ADJ  FROM  10000.  to  9384.07 

and  one  cycle  for  one  variable  has  been  completed. 

(5)  The  above  procedure  Is  repeated  for  the  second  and  all  sub¬ 
sequent  components  to  be  optimized. 

Note  that  the  same  procedure  is  repeated  for  variable  9. 


} 
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(6)  A  single  adjustment  has  now  been  made  for  each  component  for 
one  complete  search  of  the  system  component  sizes.  The  pro¬ 
cedure  is  then  repeated  for  two  more  complete  system  searches. 


(7)  The  component  whose  change  contributed  the  most  to  decreasing 
the  objective  function  is  adjusted  next  before  another  complete 
system  search  is  performed. 

(8)  The  procedure  is  terminated  when  either  no  more  improvement  in 
the  objective  function  can  be  made  (within  a  tolerance)  for  the 
component  making  the  greatest  contribution  to  decreasing  the 
objective  function,  or  the  complete  search  cycle  is  completed. 

Note  that  occasionally  no  successful  adjustment  can  be  made.  If 
the  computed  adjustment  does  not  reduce  the  objective  function,  its  value 
is  successively  reduced  to  the  original  value, testing  for  improvement  at 
a  number  of  steps  (pages  5  and  6  of  Exhibit  3). 

The  remaining  output  should  be  self-explanatory.  Remember  the  output 
is  for  two  plans  (existing  and  the  derived  system)  for  nine  flood  events 
which  results  in  18  hydrologic  simulations  at  each  control  point  and  two  eco¬ 
nomic  evaluations  at  all  damage  centers. 


SIZING  RESERVOIR  AND  PUMPING  PLANT  — 

HYDROLOGIC  PERFORMANCE  CONSTRAINED 

The  objective  for  this  example  is  to  determine  the  size  of  the 
facilities  that  will  maximize  the  system  net  benefits  while  simultaneously 
meeting  hydrologic  performance  targets  expressed  in  terms  of  desired  flow 
(storage)  target  and  corresponding  exceedence  frequency.  This  example 
extends  the  previous  example  for  the  performance  targets  of 

Exceedence  Frequency 

Reach  Target  Value  (Events  per  Year) 

1030  1200  cfs  1.0 

305  5000  ac.ft.  .05 

The  starting  values  were  selected  as  5000  ac.ft.  and  5000  cfs,  respectively. 

Pages  1  and  2  of  Exhibit  4  contain  a  listing  of  the  input  data  with 
notations  on  coding  revised  and  added.  Pages  3  through  28  contain  printout 
of  selected  pages  of  the  output. 

The  derived  optimum  sizes  are  7528  ac.ft.  for  the  reservoir  and  6044 
cfs  for  the  pumping  plant  (summary  page  28).  The  total  capital  cost  is 
$9,889,000  and  system  annual  net  benefits  are  $123,000  (benefit  cost  ratio 


of  1.15).  The  derived  values  were  adjusted  from  starting  values  of  5000 
ac.ft.  and  5000  cfs,  respectively.  The  sensitivity  of  the  solution  to 
starting  values  was  tested  by  making  a  separate  run  with  starting  values 
of  10,000  ac.ft.  and  7000  cfs,  respectively.  The  derived  sizes  were 
6,007  ac.ft.  and  6,570  cfs  costing  $9,832,000  and  resulting  in  annual  net 
benefits  of  $102,000.  The  hydrologic  performance  specified  is  achieved 
in  that  the  deqree  of  protection  provided  is  1.0  years  (protection  against 
the  annual  event)  for  reach  1030  and  .05  (protection  against  the  20-year 
event)  for  reach  305  (see  pages  15  and  26  of  Exhibit  4). 

The  output  detailing  the  progress  of  the  optimization  contains  ad¬ 
ditional  information  related  to  the  performance  target  constraints.  The 
additional  variables  are  (page  3,  Exhibit  4): 

Variable  Definition 

ISTA  =  Station  where  performance  target  specified 

INT  FLOW  =  Flow  corresponding  to  the  target  exceedence  frequency 
for  the  current  values  of  the  variables 

TRG  FLOW  =  Target  flow  for  the  target  exceedence  frequency 

FLW  OBJ  *  Component  of  penalty  applied  to  objective  function 

because  of  failure  to  meet  target  (Illustrated  later) 
for  this  station 

FLW  DEV  =  Difference  between  INT  FLW  and  TRG  FLW 

OBJ  DEV  =  Penalty  applied  to  objective  function  because  of 
failure  to  meet  target  (multiply) 

The  additional  printout  occurs  for  all  stations  where  performance 
targets  are  specified  (as  many  as  desired).  The  optimization  proceeds 
exactly  as  the  previous  (unconstrained)  example  except  that  the  objective 
function  is  penalized  whenever  the  performance  targets  are  not  met.  Note 
that  the  first  objective  function  is  extremely  large  (.951E+06)  because 
of  the  large  penalty  from  not  meeting  the  target  for  station  305  while  the 
objective  function  when  optimization  is  complete  (page  10,  Exhibit  4) 
essentially  has  no  penalty  (.106E+04).  The  computation  of  a  value  of  the 
objective  function  for  the  condition  blocked  out  on  page  5  (Exhibit  4) 
will  illustrate  the  role  of  the  penalty  assessment.  See  reference  3  for 
a  description  of  the  objective  function. 

FLW  OBJ  *  [  (FLW  DEV )/( .10  TRG  FLOW)  ]4 


Station  1030 


4 


Station  305 

FLW  OBJ  =  (— )4  =  5.988 

Objective  Function  Assessment 

OBJ  DEV  =  .0001  +  5.988  =  5.988 
0  FTN(NC)  =  (TANCST  +  ANDMG)  (OBJ  DEV  +  1 ) 

0  FTN(NC)  =  (774.217  +  265.434)  (5.988  +  1)  =  7264.80 

The  printout  at  the  bottom  of  the  pages  on  which  economic  output  is 
shown  (page  15  for  example)  summarizes  the  performance  target  and  final 
regulated  values. 

SIZING  RESERVOIR,  PUMPING  PLANT  AND  DIVERSION 


A  proposal  offered  at  past  public  meetings  has  been  to  divert  a  por¬ 
tion  of  the  runoff  from  subbasin  20  at  station  20  into  the  adjacent  water¬ 
shed  (which  is  presently  undeveloped)  both  to  reduce  flooding  in  the 
downstream  reaches  and  increase  wetlands  in  the  adjacent  watershed  to  im¬ 
prove  wildlife  habitat.  This  example  extends  the  previous  reservoir  and 
pumping  plant  example  (unconstrained)  to  include  a  diversion  from  station 
20. 

A  diversion  transfers  flow  between  locations  within  the  system. 

The  performance  characteristics  are  defined  by  a  threshold  flow  and  a 
diversion  capacity.  The  concept  of  the  diversion  relationship  is  in¬ 
dicated  in  figure  3.  Water  diverted  may  be  returned  to  the  system  at 
any  downstream  location  so  that  it  is  possible  to  characterize  facilities 
which  would  bypass  a  portion  of  flood  flows  around  a  damage  center.  Flow 
may  also  be  permanently  diverted  from  the  system,  which  will  be  done  for 
this  example.  The  cost  is  characterized  similar  to  a  pumping  plant  by 
a  capital  cost  function,  a  capital  recovery  factor  and  annual  operation, 
maintenance  and  replacement  factor. 

Table  5  (Appendix  A)  summarizes  the  performance  and  cost  data  for  the 
proposed  diversion. 

The  coding  to  include  a  diversion  at  station  20  Is  noted  on  the  listing 
of  input  data,  pages  1  and  2  of  Exhibit  5.  Note  that  it  was  necessary  to 
Include  a  dummy  reservoir  at  station  20  in  order  to  accommodate  the  require¬ 
ments  for  a  diversion. 
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Figure  3.— -Effect  of  Diversion  on  Flood  Hydrograph 

Pages  3  through  34  of  Exhibit  5  contain  selected  pages  of  the  output. 

The  derived  optimum  sizes  are  6620  ac.ft.  (index  storage)  for  the  reservoir, 
863  cfs  for  the  diversion  and  2250  cfs  for  the  pumping  plant  (summary  page 
34).  The  total  capital  cost  is  $7,099,000  and  system  net  benefits  are 
$197,000  (benefit  cost  ratio  of  1.33).  The  derived  values  were  adjusted 
from  starting  values  of  4000  ac.ft.,  500  cfs  and  1000  cfs,  respectively,  for 
the  reservoir,  diversion,  and  pumping  plant.  The  sensitivity  of  the  sol¬ 
ution  to  starting  values  was  tested  by  making  a  separate  run  with  new  star¬ 
ting  values  of  10,000  ac.ft.,  3000  cfs  and  4000  cfs,  respectively.  The 
derived  sizes  were  6648  ac.ft.,  1393  cfs  and  2160  cfs,  respectively,  costing 
$7,617,000  and  resulting  in  annual  net  benefits  of  $167,000.  In  comparison 
with  the  previously  derived  values,  it  appears  the  diversion  should  be  the 
smaller  size.  Additional  runs  demonstrate  the  value  of  testing  a  few  starting 
values  in  an  effort  to  locate  a  reasonable  optimum. 

The  hydrologic  performance  of  the  derived  system  can  be  characterized 
by  the  degree  of  protection  provided,  i.e.,  the  exceedence  frequency  of 
the  threshold  of  damaging  flows.  At  control  point  1030,  the  0  damage  ex¬ 
ceedence  frequency  was  reduced  from  about  the  five  times  per  year  event  to 
about  the  annual  event  (about  the  same  as  the  example  without  the  diversion). 
At  control  point  2030,  the  0  damage  exceedence  frequency  was  not  materially 
changed  from  the  five  times  per  year  event.  At  control  point  305,  the  fre¬ 
quency  of  significant  damage  was  reduced  from  about  the  3-year  exceedence 
interval  event  to  between  the  10  and  15-year  events.  The  residual  damages 
for  the  system  are  reduced  to  about  1/3  of  the  damages  under  existing 
conditions. 


SIZING  LOCAL  PROTECTION  PROJECTS 


Local  protection  projects  include  levees,  floodwalls  and  channel 
modifications.  Ignoring  for  the  moment  natural  valley  storage  effects, 
the  hydrologic  and  economic  effects  of  local  projects  are  truly  local. 
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i.e.,  do  not  interact  with  the  system  hydrology.  If  this  is  the  case, 
and  it  will  be  unless  the  modification  is  extensive,  then  a  local  project 
can  be  completely  characterized  performance-wise  by  a  design  Q  (or  stor¬ 
age)  and  a  flow  (or  storage)  damage  function.  Damages  are  usually  neg¬ 
ligible  below  the  design  flow  and  follow  a  curve  related  to  the  local  site 
hydraulics  and  damage  potential  above  this  point.  A  levee  or  floodwall 
essentially  truncates  the  damage  function  below  the  design  flow  (basic  hy¬ 
draulic-economic  relationship  unchanged)  while  channel  modifications  lower 
the  relationship  in  response  to  the  improved  conveyance  characteristics. 

The  concept  embodied  in  HEC-1  is  that  a  design  flow  is  associated 
with  a  unique  damage  relationship  and  therefore  if  the  range  of  feasible 
design  flows  are  known,  the  relationship  for  a  specific  design  flow  with¬ 
in  the  feasible  range  could  be  determined.  The  relationship  (flow  or 
storage-damage)  for  a  specific  design  flow  is  determined  by  interpolating 
between  the  relationships  defining  the  feasible  range.  The  relationships 
defining  the  feasible  range  are  termed  "pattern  functions;"  the  minimum 
design  damage  function  corresponding  to  the  design  flow  considered  the 
lowest  value  feasible  and  the  maximum  design  damage  function  corresponding 
to  the  design  flow  considered  the  highest  value  feasible. 

The  local  projects  considered  for  this  example  are  a  channel  modifica¬ 
tion  for  reach  1030  and  a  levee  for  reach  2030.  The  pattern  damage  func¬ 
tions  for  reach  1030  were  developed  from  water  surface  profile  and  economic 
studies.  The  minimum  design  damage  function  corresponds  to  a  "clear  and 
snag"  alternative  and  was  constructed  by  computing  water  surface  profiles 
for  a  smoothed  boundary  to  develop  a  rating  curve  at  the  index  station  that 
was  subsequently  combined  with  an  area,  elevation,  damage  relationship.  The 
design  flow  associated  with  this  function  is  1700  cfs,  the  lower  limit  of 
design  flow.  The  maximum  design  function  corresponds  to  a  40  ft.  bottom 
width,  2  to  1  side  slope  channel  enlargement  and  was  constructed  by  compu¬ 
ting  water  surface  profiles  for  modified  hydraulic  geometry  to  develop  a 
rating  curve  that  was  subsequently  combined  with  an  area,  elevation,  damage 
relationship.  The  design  flow  associated  with  this  function  is  8300  cfs, 
the  upper  limit  of  design  flow  for  the  enlarged  channel.  Table  6  (Appendix 
A)  summarizes  the  performance  and  cost  data  for  the  proposed  channel  mod¬ 
ification  for  reach  1030.  Table  6  also  contains  a  generated  damage  function 
for  a  specific  design  flow  to  illustrate  the  interpolation  concept. 

The  upper  and  lower  pattern  damage  functions  for  reach  2030  are  the 
same  and  correspond  to  existing  conditions.  The  reason  for  the  correspon¬ 
dence  is  that  the  effect  of  a  levee  is  primarily  to  truncate  the  function 
at  the  design  flow.  Some  change  is  possible  for  various  designs  if  the 
flow  area  is  greatly  restricted  by  the  levees.  The  example  assumes  no  sig¬ 
nificant  conveyance  change  from  the  levees,  though  the  methodology  does 
not  require  the  assumption.  Table  7  (Appendix  A)  summarizes  the  cost  and 
performance  data  for  the  proposed  levee  reach. 

The  existing  conditions  damage  relationships,  cost  and  runoff  hydrol¬ 
ogy  for  reaches  1030  and  2030  have  been  purposefully  made  the  same  so  that 
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the  methodology  developed  for  handling  local  projects  can  be  easily  observed. 

The  example  contains  only  local  projects  (other  damage  centers  and  alter¬ 
natives  removed)  so  that  the  difference  in  the  derived  sizes  of  the  two 
alternatives  should  only  be  due  to  differences  in  their  performance,  i.e., 
modified  damage  relationships.  A  listing  of  the  Input  data  for  this  example 
is  contained  on  pages  1  and  2  of  Exhibit  6. 

Pages  3  through  15  of  Exhibit  6  contain  sel ected  pages  of  the  output 
of  the  optimization  run. 

The  derived  optimum  sizes  are  about  5000  cfs  design  flow  for  both 
the  channel  modification  reach  and  the  levee  reach.  This  amounts  to  about 
a  0.7  exceedence  frequency  degree  of  protection.  The  total  capital  cost 
is  $207,000  and  system  annual  net  benefits  of  $30,000.  The  derived  values 
were  adjusted  from  starting  values  of  2000  cfs  design  flow  for  each  facility. 

It  is  interesting  to  note  that  while  both  facilities  began  and  ended  with  the 
same  values,  the  adjustment  route  to  the  optimum  was  different.  There  was 
no  requirement  that  they  both  end  up  the  same  size  (see  pages  3  through  5  of 
Exhibit  6).  In  addition,  note  that  while  the  values  derived  were  the  same, 
the  net  benefits  were  different  because  the  damage  relationships  were  quite 
different.  The  channel  modification  cost  $104,000  and  had  average  annual  benefits 
of  27,000  for  annual  net  benefits  of  $19,000  (benefit  cost  ratio  of  approx¬ 
imately  3.4).  The  levee  cost  $103,000  and  had  average  annual  benefits  of  $19,000 
for  annual  net  benefits  of  $11,000  (benefit  cost  ratio  of  approximately  2.4). 


SIZING  RESERVOIR,  PUMPING  PLANT,  DIVERSION, 
AND  UNIFORM  PROTECTION  LOCAL  PROJECTS 


This  final  example  includes  all  the  proposed  components  that  have  been 
previously  illustrated.  The  optimization  will  be  unconstrained  and  the 
uniform  protection  level  option  for  the  local  projects  will  be  used.  The 
uniform  protection  level  option  will  in  effect  cause  a  "degree  of  protec¬ 
tion"  to  be  optimized  for  the  two  local  protection  projects.  A  complete 
listing  of  the  input  data  is  contained  on  pages  1  through  3  of  Exhibit  7 
and  the  complete  output  on  pages  4  through  39. 

The  derived  optimum  sizes  are  6701  ac.ft.  for  the  reservoir,  0.2  ex¬ 
ceedence  frequency  for  the  levee  and  channel  projects  (2947  cfs  for  the 
channel  modification  and  7660  cfs  for  the  levee),  670  cfs  for  the  diversion 
and  2450  cfs  for  the  pumping  plant  for  a  total  capital  cost  of  $7,408,000 
and  system  net  benefits  of  $196,000  (benefit  cost  ratio  of  1.31).  The  op¬ 
timum  sizes  were  adjusted  from  starting  values  of  4000  ac.ft.  for  the 
reservoir,  0.2  exceedence  frequency  (uniform  protection)  for  local  projects, 
500  cfs  for  the  diversion  and  1000  cfs  for  the  pumping  plant.  A  comparison 
of  Exhibits  5  and  7  Indicates  that  the  inclusion  of  local  projects  has 
very  little  effect  on  the  optimum  sizes  of  the  major  facilities  (reservoir 
and  pumping  plant).  The  diversion  capacity  was  lowered  slightly  from  that 
derived  in  Exhibit  5  which  probably  means  that  it  is  more  efficient  to  pro¬ 
tect  reach  2030  by  the  levee  project. 
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OBJECTIVE  OF  THE  FLOOD  CONTROL  SYSTEM 
COMPONENT  OPTIMIZATION  UTILIZING  HEC-1 


The  optimization  algorithm  (or  search  procedure)  discussed  In  this 
training  document  has  been  developed  to  assist  the  planner  In  systematically 
and  efficiently  screening  a  large  number  of  possible  flood  control  alter¬ 
natives.  Although  there  Is  an  upper  limit  to  the  number  which  can  be  satis¬ 
factorily  and  economically  optimized  in  one  particular  computer  run.  It  Is 
still  possible  to  analyze  a  large  number  of  components  by  grouping.  In  the 
Phoenix  Urban  Study,  Los  Angeles  District  Corps  of  Engineers  (reference  6), 
there  were  eight  upstream  storage  alternatives  to  be  evaluated.  Although 
each  component  was  analyzed  Individually,  It  was  possible  to  determine  which 
component  and  combination  of  components  were  economically  feasible  by  making 
several  runs  In  groups  of  two  and  three  components  and  comparing  the  economic 
and  hydrologic  consequences. 

It  should  be  emphasized  that  the  optimization  procedure  of  HEC-1  Is 
a  planning  tool  for  determining  potential  and  economically  feasible  flood 
control  alternatives.  Once  those  that  have  potential  are  selected,  then 
a  more  detailed  simulation  of  the  operational  and  hydraulic  characteris¬ 
tics  of  a  particular  component  will  probably  be  required  as  various  stages 
of  study  (leading  to  design)  are  undertaken. 
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TABLE  1 

HYDROLOGIC  DATA 
(Existing  Conditions) 


DRAINAGE  AREA 


Area 

Subbasin  (square  miles) 


10 

35.1 

20 

35.1 

30 

10.0 

TOTAL  80.2 

SUBBASIN  RUNOFF 

SYNTHETIC  STORM  EVENT 

(hourly  values) 

Inflow  to 

Inflow  to 

Inflow  to 

Sta. 

JOJcfs) 

Sta. 

20  (cfs) 

Sta. 

30  (cfs) 

24 

2200 

24 

2200 

8 

730 

24 

1840 

24 

1840 

8 

615 

26 

1540 

26 

1540 

9 

515 

33 

1250 

33 

1250 

11 

415 

50 

995 

50 

995 

17 

330 

85 

775 

85 

775 

28 

255 

190 

605 

190 

605 

63 

200 

375 

470 

375 

470 

125 

155 

515 

365 

515 

365 

170 

120 

590 

280 

590 

280 

195 

93 

660 

215 

660 

215 

220 

72 

710 

160 

710 

160 

230 

54 

760 

120 

760 

120 

255 

41 

800 

95 

800 

95 

265 

32 

840 

77 

840 

77 

280 

26 

910 

66 

910 

66 

305 

22 

1040 

59 

1040 

59 

350 

20 

1290 

53 

1290 

53 

430 

18 

1920 

49 

1920 

49 

640 

16 

3000 

42 

3000 

42 

1000 

14 

3950 

40 

3950 

40 

1320 

13 

4600 

38 

4600 

38 

1540 

12 

5080 

35 

5080 

35 

1650 

I  h 

II 

5360 

33 

5360 

33 

1800 

I  V 

II 

5370 

30 

5370 

30 

1810 

1  1 

11 

5100 

30 

5100 

30 

1690 

10 

4600 

29 

4600 

29 

1530 

10 

3980 

27 

3980 

27 

1330 

g 

3330 

25 

3330 

25 

1110 

* 

a 

2720 

25 

2720 

25 

900 

* 
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TABLE  1  (Continued) 


HYDROLOGIC  DATA 
(Existing  Conditions) 


Reach  10-30  Hod.  Puls  Routing  Criteria1 


Storage  (ac.ft.) 

0 

SO 

475 

940 

2135 

3080 

6300 

Outflow  (cfs) 

0 

200 

1020 

2050 

6100 

10250 

24000 

Reach  20-30  Mod  Puls  Routine  Criteria1 
Storage  (ac.ft.)  0  50  475 

940 

2135 

3080  6300 

Outflow  (cfs)  0  200  1020 

2050 

6100 

10250  24000 

Outflow  Culvert  (Sta 

Storage  (ac.ft.) 
Outflow  (cfs) 


30)  Hod.  Puls  Routing  Criteria1 

0  400  1000002 

0  1200  1200 


1/ 


Storage-outflow  data  should  extend  beyond  the  maximum  values  computed  In 
the  multi  flood-multi  pi an  options. 


Zj 

Note  that  the  outflow  becomes  constant  and  equal  to  1200  cubic  feet  per 
second  when  the  detention  storage  equals  or  exceeds  400  acre  feet. 


TABLE  2 

ECONOMIC  DAMAGE-FREQUENCY  DATA 
(Existing  Conditions) 


Exceedence  Frequency 
(Events  per  Yr) 

Damage  Center  1030 

Type  1 

Flow  Damage 

(cfs)  ($1000) 

Type  2 
Damage 
($1000) 

Type  3 
Damage 
($1000) 

1 

6.000 

1030 

0.00 

0.00 

0.00 

, 

5.500 

1130 

0.00 

0.00 

0.00 

4.500 

1380 

0.10 

0.50 

1.00 

;  ' 

* 

3.500 

1740 

0.20 

0.70 

1.50 

| 

i 

2.500 

2280 

0.30 

1.50 

3.20 

T. 

*  t 

1.500 

3200 

0.30 

2.20 

4.70 

s 

.900 

4220 

0.40 

2.90 

6.50 

*  1 

•  i 

.700 

4800 

0.50 

3.50 

7.80 

i  i 

1  ] 

.500 

5620 

0.60 

4.00 

9.30 

i 

.350 

6480 

0.70 

4.70 

11.00 

.250 

7340 

0.80 

5.80 

13.70 

i 

.150 

8540 

0.90 

6.60 

15.60 

.100 

10000 

1.00 

8.00 

19.00 

1 

.050 

12100 

1.20 

10.30 

23.00 

.  i 

.020 

15100 

1.50 

15.00 

27.80 

.005 

21000 

1.80 

18.10 

30.20 

i  | 

i 

Damage  Center  2030 

! 

j 

Type  1 

i 

i 

Exceedence  Frequency 

Flow 

Damage 

(Events  per  Yr) 

_ tiisl _ 

telQQQ) 

6.000 

1030 

0.00 

5.500 

1130 

0.00 

4.500 

1380 

1.60 

i 

3.500 

1740 

2.40 

| 

2.500 

2280 

5.00 

1.500 

3200 

7.20 

i 

.900 

4220 

9.80 

.700 

4800 

11.80 

I 

.500 

5620 

13.90 

1 

.350 

6480 

16.40 

i 

i 

.250 

7340 

20.30 

1 

.150 

8540 

23.10 

i 

i 

.100 

10000 

28.00 

1 

.050 

12100 

34.50 

1 

l 

.020 

15100 

44.30 

.005 

21000 

50.10 

t 


TABLE  2  (Continued) 
ECONOMIC  DAMAGE-FREQUENCY  DATA 
(Existing  Conditions) 


Exceedence  Frequency 
(Events  per  yr) 

Damaae  Center  305^ 

Storage 

(ac-ft) 

Type  1 
Damage 
($1000) 

Type  2 
Damage 
($1000) 

.700 

1500 

0.00 

0.00 

.600 

2300 

37.50 

10.50 

.450 

4000 

75.00 

15.00 

.250 

7000 

1125.00 

52.50 

.100 

12500 

3150.00 

105.00 

.050 

20000 

5850.00 

202.50 

.020 

28000 

7050.00 

300.00 

.010 

37000 

9000.00 

390.00 

.005 

50000 

10650.00 

540.00 

.002 

76000 

11250.00 

585.00 

Flood  Ratios  for  Multiflood,  Multiplan  Evaluation 
0.25  0.30  0.50  0.70  1.00  1.50  2.20  3.25  4.40 


Note  that  the  damage-frequency  relationship  (for  damage  center  305)  Is  a 
function  of  storage  and  not  discharge. 


TABLE  3 

RESERVOIR  PERFORMANCE  AND  COST  DATA 


Low  Level  Outlet 


Area  of  Opening 
Orfice  Coefficient,  C, 

In  the  general  expression 

Q=  C  A  (2gH)Exp* 

(free  discharge) 

Cenderllne  Elevation  of  Orfice 
No  Tallwater  (no  submergence) 
Exponent  of  head  (Fxp.) 


Overflow  Spillway 

Type 

Length 

Weir  Coefficient,  C, 

in  the  general  expression 

Q=  C  L  H3/2 


Cost  and  Site  Characteristics^ 

Capacity  (ac.ft.)  0  2500  4000 
Elevation  (ft)  965  1000  1015 
Cost  ($1000)  0  1500  2400 


=  35  ft2 


*  0.71 

*  975  ft 

■  0.5 


=  Ogee 
*  35  ft 


*  2.86 


5200 

6800 

9000 

11500 

15500 

21000 

1030 

1045 

1060 

1075 

1090 

1105 

3000 

3600 

4350 

4950 

5550 

6000 

Annual  Cost  Data 


Annual  Operation  and  Maintenance  =  2.33  of  Capital  Cost 
Discount  Factor  (Capital  Recovery)*  5.043 


Constraints 


Reservoir  size  must  be  in  range  of  0  to  25,000  ac.ft. 


30000 

1120 

7200 


1/ 


Capacity-elevation  data  should  extend  beyond  the  maximum  values  computed 
in  the  multi  flood-multi  pi  an  options  and  the  maximum  reservoir  size  designated 


TABLE  4 

PUMPING  PLANT  PERFORMANCE  AND  COST  DATA 


V. 


u. 


Cost  and  Performance  Data 


Capacity  (cfs)  0 

250 

500 

1000 

2000 

6000 

8000 

10000 

Cost  ($1000)  0 

670 

1000 

1600 

2300 

6000 

7860 

8670 

Annual  Cost  Data 

Annual  Operation  and  Maintenance  =  2.3 %  of  Capital  Cost 
Discount  Factor  (Capital  Recovery)  s  5.04%  j 
Annual  Power  Cost  =  $100,000 


Sizing  and  Operation  Data 

Pumping  plant  must  be  between  0  and  10,000  cfs. 

Pumps  activate  at  storage  level  (at  station  30)  *  1500  ac.ft. 


y 

Annual  power  cost  is  adjusted  based  on  the  difference  In  computed 
volumes  at  the  pumping  facility  as  system  component  sizes  vary  from 
specified  initial  values  to  optimized  values 


TABLE  5 

DIVERSION  PERFORMANCE  AND  COST  DATA 


Performance  and  Cost  Data 


Capacity  (cfs)  0 

1250 

2500 

3750 

5000 

7500 

10000 

15000 

Cost  ($1000)  0 

1500 

2600 

3400 

4200 

5200 

6100 

7500 

Annual  Cost  Data 

Annual  Operation  and  Maintenance  3  1.5%  of  Capital  Cost 
Discount  Factor  (Capital  Recovery)  *  5.04% 


Operation  and  Constraints 

Diversion  activation  threshold  *  1,500  cfs 

Size  limit  between  0  and  20,000  cfs 


A-7 


20000 

8300 


TABLE  6 

CHANNEL  MODIFICATION  COST  AND  PERFORMANCE  DATA 


A-F 


V  _  * 


TABLE  6  (Continued) 

CHANNEL  MODIFICATION  COST  AND  PERFORMANCE  DATA 


Performance  and  Cost  Data 

Capacity  (cfs)  1700  5000  5500  7000  8300  9300 

Cost  ($1000)  42  103  149  222  283  340 


Annual  Cost  Data 

Annual  Operation  and  Maintenance  *  2.3%  of  Capital  Cost 
Discount  Factor  (Capital  Recovery)  •  5.04  % 


Design  Limits 

Minimum  Design  Q  *  1700  cfs 
Maximum  Design  Q  ■  8300  cfs 
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TABLE  7 


LEVEE  COST  AND  PERFORMANCE  DATA 
Damage  Center  2030 


Flow 

(cfs) 

Minimum  Design 

Damage  Function 

Damage 

($1000) 

Maximum  Design 
Damage  Function 
Damage 
($1000) 

1030 

0.00 

0.00 

1130 

0.00 

0.00 

1380 

1.60 

1.60 

1740 

2.40 

2.40 

2280 

5.00 

5.00 

3200 

7.20 

7.20 

4220 

9.80 

9.80 

4800 

11.80 

11.80 

5620 

13.90 

13.90 

6480 

16.40 

16.40 

7340 

20.30 

20.30 

8540 

23.10 

23.10 

10000 

28.00 

28.00 

12100 

34.50 

34.50 

15100 

44.30 

44.30 

21000 

50.10 

50.10 

Performance  and  Cost  Data 

Capacity  (cfs) 

1700 

5000 

5500 

7000 

8300 

9300 

Cost  ($1000) 

42 

103 

149 

222 

283 

340 

Annual  Cost  Data 

Annual  Operation  and  Maintenance  =  2.3%  of  Capital  Cost 

Discount  Factor  (Captial  Recovery)  =  5.04% 


Design  Limits 

Minimum  design  Q  =  1700  cfs  Maximum  design  Q  =  8300  cfs 


EXHIBIT  1 
HYDROLOGIC  MODEL 
(Existing  Conditions) 


'LOUD  C'INT»UL  System  COMPONENT  OPTIMIZATION 

»»n»ntnf.ic  MnoEt 
Existing  conditions 
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GRAVITY  outlet  through  levee 


EXHIBIT  2 


MULTIFLOOD,  MULTIPLAN  MODEL 
(Economic  Evaluation  of  Existing  Conditions) 
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Pl»«  FLOW  »*D  STOWAGE  (END  OF  PERtOO)  SUMMARY  FOR  MULTIPLE  PL»N-R»TIU  ECONOMIC  COMPUTATIONS 
FLU«S  IN  C  US  I C  FEET  PER  SECOND  (CUPIC  MEIERS  PER  SECOND) 

AREA  IN  SQUARE  MILES  (SSUARE  MILOMETERS) 
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